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SUMMARY 


A method is described for the solution of the nonlinear equations 
for supersonic conical flow. The procedure is mostly a numerical one 
based on the method of characteristics and the relaxation process. A 
procedure for calculating the position of the shock is inherent in the 
analysis. The method is applicable to any conical flov. 

As an illustration, the flow about a triangular wing with super- 
sonic edges is presented. 


INTRODUCTION 


A knowledge of the characteristics of bodies in supersonic flight 
is depeMent on both theoretical and experimental information. Most 
theoretical aerodynamic data are obtained by the use of linearized 
theory which permits only very slender bodies and small angles of 
attack. In determining the value of such information and the practical 
range over which the parameters such as body diameter, wing thickness, 
and angle of attack may be allowed to vary, one must rely on con 5 )arison 
of these results with either experiment or a nonlinear theory. Because 
of the limited experimental facilities available and the expense 
involved in their operation and in model construction, the existence of 
a practical nonlinear theory becomes iu 5 )ortant. Further, if more accu- 
rate results are desired than may- be obtained from a linear theory, a 
nonlinear one is required. 

In the case of conical flows, the solutions which have been obtained 
up to the present time have been found mostly by the use of linearized 
theory. Busemann (reference l) introduced this procedure. References 2 
and Z} "to mention only two Investigations, present extensive studies of 
thin conical wings. Browne, Friedman, and Hodes (reference 4 ) have 
solved a special wing-body problem. In reference 5^ ^ fairly general 
method of solution is given for fuselage-type conical bodies. Linear- 
ized solutions, satisfying the exact boundary conditions, have been dis- 
cussed by Lapoidie and Bartels (reference 6). Moore (reference 7 ) and 
Broderick (reference 8) have obtained second-order linearized solutions 
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for several conical flows. More general analyses applicable to conical 
problems include, for exan^jle, Eward’s work given in reference 9 (this 
is applicable in the lifting case only wben there is a supersonic 
leading edge) and Spreiter’s solutions (reference 10 ) for very slender 
bodies and low-aspect- ratio wings. 

The application of nonlinear theory to conical flows has been lim- 
ited mostly to Taylor and Maccoll’s solutions for the ri^t circular 
cone at zero angle of attack (references 11 and 12 ). Solutions for the 
yawed circular cone are given in references 13 to 15. Busemann (refer- 
ence 16) and Ferrari (reference 17) have presented a general discussion 
of the nonlinear problem. Ferri (reference I8) has given some general 
analysis of the problem with particular reference to the yawed cone. 
However, until the present, no attemrpt has been made to present a 
detailed method for the solution of the general conical- flow problem. 

The solution presented herein is of necessity mostly a numerical one 
using the methods of characteristics and relaxation. In the latter 
case, reference to work by Southwell and Emmons (see, e.g. , refer- 
erences I9 and 20 ) has been found valuable. A discussion, by Li^thill, 
of the shape of the shock about a conical obstacle (reference 21 ) also 
proved useful. 

This investigation was conducted at Brown University under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 

The author wishes to express his appreciation to Dr. W. D. Hayes 
for many valuable suggestions received in the course of discussions with 
him and to Dr. G. F. Carrier who read the manuscript critically and 
carefully. 


SYMBOLS 


A rotation, defined by equation (lO) 

c velocity of sound divided by maximum velocity 

c* velocity of sound 

Cp pressure coefficient, difference between local and free- 

streaon static pressures divided by free-stream dynamic 
pressure 

J • mechanical eq.ulvalent of heat 

function of velocity and its derivatives, defined by 
equation (Uo) 


K 
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M 

%1 

%2 

Q 

R, e, cp 
R 

S = JS‘/ 7 R 

S* 

%1 

%2 

T 


u, V, V 


u 


V’ 

V 


a 

P 


Mach number of free stream 

con5)onent of free-stream Mach number taken normal to vlng 
leading edge 

con^jonent of free-stream Mach 'number taken behind and 
normal to shock 

function of shock position, slope, and free-stresim Mach 
number, defined by equation ( 33 h) 

spherical cooidinates 

gas constant 

entropy 

speed normal to wing leading edge 
speed behind and normal to shock 
speed tangent to wing leading edge 

velocity conqponent lying on surface of sphere throu^ 
conical vertex (\/v 2 + w 2 ^ 

ratio of velocity con^ionents in'r-, 0 -, and cp-dlrectlons, 
respectively, to maximum velocity 

ratio of free-stream velocity to maximum velocity 

ratio of Cartesian velocity con5)onent8 to speed of sound 

velocity vector 

V’ divided by maximum velocity 
maximum velocity 

ratio of velocities in and cp-directions (-v/w) 
streamlines ■ 


7 


adiabatic exponent 
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e 

V 

Subscripts: 
9, cp 


arc--tangent of slope of shock , on surface r = Constant 
slope of hody-tangent plane, defined in figure k 
angle defined in figure k and equations {2k) 
eq.uatlon of body- 

shock angle, measured in plane lying normal to leading 
edge of surface to -which shock is at-tached 

angle defined in figure k and equations {2k) 

angle defined in figure ^ and equations {2h) 

differentiation -vriLth respect to 0 and qp, respectively 


AKALYSIS 


The equations of motion of a perfect isoenergetic fluid may be 
-written according to reference 22 as 


and 


-where 



( 1 ) 


( 2 ) 


(c’)^ = - V« . V«) (3) 

and -where V* is the velocity vector; S*, the entropy; J, the mechein- 
Ical equl-ralent of heat; 7 , the adiabatic exponent; R, the gas con- 
stant; the ma-xi mum speed obtainable on eacpanslon into a vacuum 

(also constant); and ■ c*, the velocity of sound. A change of -mrlables 
may now be made to nondlmensionalize the quantities. Let 
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s = JS‘/rR 
V = V‘/¥max r 
c = c*/ Wmay 


( 1 ^) 


Using these new variables, equations (l) and (2) were expanded in spher- 
ical coordinates r, B, and cp (fig. l). This was done for conical 
flow (so that neither the velocity' nor entropy are functions of r). 
Noting that u, v, and w are the conponents of V in the r-, 0-, 
and cp-directions, respectively, one obtains 


/ 

u(v - Ug) - w(vg - + » cot ej = C% 


■(' 


V Wa - 


'9 


cot e) - JJ!sl - v) = 

/ \sin 0 / sin 0 


0 sin 0 


+ w 


and, using equation (5a-), 


sin 0 + (v sin 0)„ + w^ = 
sin 0j_ e 9j 


u(v2 + w?) + |(t 2 + + _^_(v2 + 


(5a) 

(5b) 


(5c) 


( 6 ) 


Equations (5) yield 


VSg 


+ ^ 
sin 0 


= 0 


(7) 


This says sliiq)ly that the entropy is constant on the intersection 
of a stream surface with a sphere through the vertex of the conical 
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flow. Henceforth, such an Intersection will he termed, for convenience, 
a streamline. Denote hy p = p(9) the equation of such a streamline. 
Its equation, in terms of the velocity components, is 


^0 + 


sin 0 


= 0 


(8) 


Thus the entropy is a function of p only and will hence he determined 
solely hy conditions at the juncture between the disturbed and 
undisturbed flow. 

. Solving equations (5a) and (8) for v and w as functions of p 
and u, one obtains 


■where 


and 


V = U0 - A 


= -J^-oA 

sin 6 


> 


(9) 


a 


V _ ^qp 

w " Pg sin e 




A = 




0 




sin 0 


1 + 


( 10 ) 


The A is essentially the rotation. Then, using these relations, the 
•two remaining ones (equations (5b) and C6)) become 


uA + 




(oA sin 9 )q 

\sin 0 / 

sin 0 

sin 0 



( 11 ) 
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sin^0 


) u,,fi . 

I + Wp I 

-J^S 

ro 

/ c2 / 

f J Prt ' 

sin'^0 

^ c^ sin%y 


Uq cot 9 




bIu^gJ sin^e 


—i— f^A sin 9) + (oA)”] + 

sin 0 [_ '9 ^ 


-2UgUgg + 


(ii?e ® - 2\e) - Ve“e ' 




. -u(l + a2) + uUq + ( + auja- + /2u 

_ \sin 0 0/0 sin 0V ^ 


u cot 0 

fp 


|5^(a>,>“<p + “Vp) - + 


sin 0> 


Equation (12) can te written in this alternative foim: 

1 _ fa - A)" _ (i^ - ■^)^~ | , , _ r _ (u, - Af 

L c2 c2 J ' ® ^2j 


\sin 0 4 


sin 0 




isin 0 ) 


- a) cot 01+ 


ksin 0 


c^ sin 0 sin 
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It is noted that. If the flow is irrotational, then Sp = 0, 

A = 0, and eq.uation (U) is identically satisfied, and the ri^t-hand 
side of equation (12) is also 0. In such a case, one need merely find u 
satisfying the left-hand side of equation (12) set equal to 0 and, of 
course, obeying the boundary conditions. It is worth noting that • u 
jnay, in this case, he thou^t of as a velocity potential (actually, the 
velocity potential would he ru(0,cp)). 

To solve a conical-flow problem eq.uatlons (U) and (12) or (I 3 ) 
need to be solved, subject to the appropriate boundary conditions. There 
are two main difficulties which will arise. One comes from the boundary 
conditions from a shock surface, the location of which is initially 
unknown. The other difficulty is caused by the faci that, in general, 
the differential equations form a set which in every case has elliptic 
regions and in many cases will also have hyperbolic ones. The condition 
for the system to be elliptic is that 

c2 > v2 + (l4) 

One may note that in the case of small disturbances, such as linearized 
flow, this condition corresponds to ■whether the point in question is 
inside or outside the free-stream Mach cone extending do-wnstream from 
the vertex:. In configurations such as a body of revolution at angle of 
attack, the domain inside the shock -will be entirely elliptic, while in 
the case, for example, of a wing -with supersonic leading edges, as will 
be discussed in more de'tail later, there -will be both elliptic and 
hyperbolic regions. Any part of a hyperbolic region which is con5)letely 
Independent of the elliptic region may be conveniently calculated by the 
method of characteristics, while any other region, either elliptic or 
hyperbolic, will be found in this report by a relaxation process. 

To set up the method of characteristics for the systems, it is 
noted that there are two second-order quasi- linear partial differential 
equations, equations (U) and (13)^ for p , and . u, -where their second 

■ derivatives appear only as the first derivatives of a. 

Note that these two quantities are 

a = 


T = 
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Then eq.uations (U) and ( 13 ) may he rewritten as 




- ■" ^ ‘=°'^ ® 


') 


Eind 


yjl + g^ /v^ _ y _ ^cp \ _ 

“■ \c^ 1 + a^/ \ ® cc sin 6 / 


V cot 9 


The auxiliary equations required are 

dT = Tg de + T(p dtp 
da = ttg d0 + Ojp dcp 

The characteristics of this system are given hy 


sin 0 dqp = 


d0 


1 . v£V c2 


^ ± , /i£ + _ 1 


2 2 

c c 


Also^ the streamlines are given hy 


d0 + a sin 0 dcp = 0 


( 16 ) 


( 17 ) 


( 18 ) 


( 19 ) 


( 20 ) 


From equations (16) throu^ (18), the condition to he fulfilled hy da 
and dT on the characteristics is also obtained. This is, for the 
characteristics defined hy equation (19), 
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The ± signs in equations (I 9 ) and (21) coincide. Thus, referring to 
figure 2 , to continue a solution from a curve ahc (not a characteristic 
curve) to a point d, one constructs the straight segments ad and cd 

hy equation ( 19 ) and hd hy equation (20). Then equation (21) enables 
one to find a and T (and hence v and w hy equation (15)) at 
point d. Then u is obtained by 


du = Uq d© + u^ dq) 

or, along the Streamline bd, using equations (15) and (20), 

du = ^T sin 9 \Jl. + a^)dcp 

Finally Pg is needed. Since bd is a streamline, = P^^. 
the characteristic ad. 


( 22 ) 


Then, on 


dP = Pb - Pa 

= Pg d© + P^p (hp 
= p„(d© + a sin © dcp) 

C7 




(23) 
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■wiilcli gives at d and this con 5 >letes the determination of the con- 

ditions at point d. Recall that, of course, Sp is constant along hd. 

In general, the construction of the solution in the hyperbolic 
region hy the method of characteristics is complicated hy the fact that 
the houndaiy from -which the solution is begun will be an initially 
unknown shock surface. Construction of the shock solution may be carried 
out in a manner similar to that given, for exanqple, in reference 23 , for 
the case of plane or axially symmetric flow. 

Referring to figure 3, suppose "that the shock is to start at 
point D. Its initial angle may be found by assuming that near D there 
is a plane -two-dimensional shock. This may be easily calculated from 
the usual oblique- shock theory (see, e.g., reference 24), idiere one 
must first calculate the con 5 )onents of the free-stream flow normal and 
parallel to the ray passing from -the vertex of the cone throu^ the 
point D. Then the conqionent of -the Mach number normal to this ray is 
used as the initial Mach number in the oblique-shock relations with the 
angle of flow deflection also measured peipendicular to that ray. One 
can then determine conditions immediately behind the shock. So long as 
the shock remains plane, the Cartesian veloc:J.ty conqjonents, in the 
region influenced only by the shock, -will remain constant. 

Thus, referring to figure 4(a) suppose that the boundary surface to 
•which the shock is attached is tangent to QAH at its leading edge QA. 

The shock will be tangent to the plane QAJ at QA. One then seeks to 
determine the equation of the shock tangent to plane QAJ and the 
Cartesian velocity conq)onents immediately behind the shock, in terms of 
the free-stream Mach number M, the position of the leading edge 
and the slope of the body- tangent plane g . Thus, 


Ml 


R1 


= M ^1 - cos cos^ 


(24) 


Tj = ZHKC = tan"^ 


sin p. + tan -v tan e 


cos (X cos -v(tan v - sin p. tan e) 


where is the con^onent of the free-stream Mach number normal to 

the leading edge and q is the angle of flow deflection throu^ the 
shock. Using 4 ol’tain from oblique- shock tables 

(e.g., reference 24) the shock* angle 0^^ (angle JKC in fig. 4(a)) and 
■the Mach number li ^2 behind the shock, normal to it. Then the velocity 

component normal to the shock, behind it, is 
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SN2=\/(l-Sr^) 


Sm = cos COB V 


1 + 


' 2 

1 +2Ljl1m2 
2 


is the tangential con^onent. 

Then the Cartesian velocity can5>onents Vy, and behind 

the shock are 


Vy > = Sqi ■< sin V 


cos p. cos V 


cos V- sin [JL 


^2 ^ 

I 2 ^ ^ 

/I - cos H cos V 


1 - cos ^ cos^ 

-cos v(sin n "tan t] + cos n sin v) 
sin V tan t] - sin [j. cob (x cob^ 


The eq,uation of the plane portion of the shock is, in spherical 
coordinates. 


tan 0, 


tan 0 = 


I - cos^ cos^V 
cos li. cos V 


^tan 0^ sin p. cos v + ^V os cp - (tan 0„ sin v - tan p)sln cp 
^ ” cos p/ ' ” 
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Now consider again figure 3« Suppose that the shock is known up 
to a point A and the conditions at a point B, lying on a characteristic 
from A, are also known. One then seeks the conditions at point C on 
the shock. The position of C may he found from the intersection of the 
characteristic from B, found from equation (I9), and the line segment AC, 
having the slope tan 8^. Then, the conditions at C may he determined 

from the differential expression (21) along BC and, the conditions at the 
shock. Let subscripts 1 and 2 refer to conditions upstream (to the 
left of the shock in fig. 3) and downstream of the shock, respectively. 
Then, since the conponents of velocity tangent to the shock are contin- 
uous across it, there hold 


2 sin 8 - w^ cos 8^ = ^v^ sin 8 - w^^ cos 8 

and, for the Junp in normal velocity across the shock (recall that u, 
V, and w are the ratios of physical velocity to ma xim um velocity), 
there holds 



^v^ cos 8 + w^ sin cos 8 + w^ sin 8) 


7-1 
7 + 1 


1 - - ^v^ sin 8 - w^ cos 5^^ 


(29) 


But, because u^^, v^, and w^^ are siuply the ' conponents of the free 

stream which lies parallel to the 0 = 0 axis, one obtains. 


where 


Uj^ = U COB 0 

v^ = -u sin 0 





(30) 




Ik 


NACA IN 2651 


Then, the solution of equations (28) to (30) yields 


Ug = u cos 0 


02 = 


•vdiere 


= -cot 8 j 


' Q + sin S ' 
iQ - COS%; 




= u^ sin^e/ ' ^ + sin^'^ 

vcos'^ / 


Q =. 


/I + Z-ZJl sin^e cos^S ’ 
7-1 2 

7 til T^M^ein^a 


hence, one may consider that, from equations (32) and (33) ^ 




or estpanding in a Taylor’s series along the characteristic BC (fig. 


(|)j(°c - °b) + 


(31) 

(32) 

(33a) 

( 33 h) 

( 3 k) 

3 ), 


Tc = Tb + 


( 35 ) 



MCA TN 2651 


15 


Then, using equations (21) and (35), there is obtained 



Hence one finds from equation ( 36 ), then T^ from equation (35), 

and hence v and w. The quantity 8 q follows from equation (32), 
u^ follows from equation ( 31 )^ and the remaining quantities are 

obtained as in the general case already discussed. This procedure wi2JL 
enable one to calculate the shape of the shock and the physical quan- 
tities on it. However, the reader is cautioned to note that in most 
cases the shock will be in nearly the same direction as the character- 
istics imme diately downstream' of it. Hence the point B should be taken 
close to point A. 

In some cases, such as the low-pressure side of a lifting wing, 
there will be no shock, but rather an expansion of the flow. In such 
a case the solution is -sir^jlified because the flow is irrotational nnfl 
the velocities vary continuously. The solution may be started by a 
Prandtl-Meyer expension in a manner analogous to that used where there 
was a shock. As in the shock case, the nonual conq)onent of the Mach 
number Mj^^ and the angle of flow deflection t] are found by equa- 
tions (24). Then,- referring to figure 4(b), angle CKJ^ is the Mach 
single corresponding to Hence, 


^ CKJ^^ = sln-l(l/t%J 


(37) 
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This defines the plane of the start of the expansion fan (OKAJ^^^ . Then, 
using as the initial Mach numher and t) as the deflection angle, 

one can immediately obtain from tables of Prandtl-Meyer flows (see 
reference 2k) the normal Mach numher 1^2 "behind the flow. Then 

and Sjj2 are found by equations (25) and Jhe Cartesian velocity com- 
ponents, by eq.uation (26). Then single tj + CKJ2 is the Mach angle 
corresponding to 1%2" Hence 


^ CKJg = sin“l(l/%2) - T (38) 

The eq.uations of the planes of the start (OKAJ^^) and end (0KAJ2) of the 

fan of the expansion are then found by using equation (27), where, 
instead of 0^, one uses angle CKJ]^. (equation 37) and angle CKJ2 (equa- 
tion (38)) for the beginning and end, respectively. To find the veloc- 
ities inside the fan, sinq)ly take a series of deflection singles between 
0 and q and calculate, as already described, the corresponding v^^, 

Vy, and v^, and the angle, corresponding to angle CKJ2, of the plane 
where those velocities apply. 

Since the end of the expansion fan is not a characteristic of the 
flow, the solution may bq continued ;from it in the manner already 
described. . 

The remaining kind. of. point which can arise is at the surface of 
the body in the flow (in fig. point E). Here there is specified the 
equation of the body, ,0g = 0g(cp). Hence 



and therefore, in equation (21), da is known and dT may be solved 
for Immediately. This conpletes the solution of a purely hyperbolic 
region. 

If the flow in a given domain is such that the system (equa- 
tions (11) and (12)) of differential equations defining the flow is 
elliptic, the solution must be found in some other manner. Also, if a 
hyperbolic region is bounded by an unknown elliptic region, that hyper- 
bolic region will interact with the elliptic one so that both must be 
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solved concurrently. In either case, the relaxation procedure (see, e.g., 
references I9 and 20 ) appears to he the hest for this portion. In most 
conical- flow problems, the shock will he fairly weak, at least as com^ 
pared with a normal shock, because, it must remain attached to the body 
in order to have conical flow. Furthermore, even if it is of appre- 
ciable strength, that strength will probably vary fairly slowly about 
the body. Hence, except at singular points (see reference I8), it 
appears that the entropy gradients will be small and, to a first approxi- 
mation, may be ignored. In such a case, one need solve only equa- 
tion (12) with the rlgh’t-hand side set equal to 0 . The effect of rota- 
tion may then be found by an Iteration process. For this, u, u^j, a, 

Sp, and in equation (U) are as found for the irrotational case 

and A is then determined. These values determ in e then the ri^t-hand 
side of equation (12) and the left-hand side may then be relaxed again. 
The process may be repeated if necessary. 

To apply the relaxation procedure to the' system, one must first 
write the differential equation in finite difference form. Referring to 
figure 5 ^ at point E there holds approximately. 


\ ®)e 2A0 


(-^) = 
\sin eJ -6 


~ ^F 
2^ sin 


(“ee)® = 


~ ^E ^ 


> 


( ) _ ^A - 

(sin 0/-g ^ bIix 


cpcp 


vsln (Zitp sin 


m 


However, in solving the (irrotational) problem, it appears easiest for 
conqjuting purposes to relax u by using u^^ and u^^ in finite dif- 
ference form and to correct periodically the other terms. Thus one gets, 
for point E 
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(^3 - 2i^ + 



+ ( 1 ^ - 2ug + up)| 


(c sin 0 )" 

r 

sin 0 )^ 


= -K (4l) 


■where 


K = u 2 - 


u. 


ucp' 


+ Uq cot 01+ 




^ ''V'^ ‘ c2^20y 


- 2 


.^0^(pV 


sin^0 


Also it is seen that, since the difference equation (4l) has variable 
coefficients, there is no particulax ad-vantage in mfi'k-l ng A0 and ^ 
everywhere the same. In particular, the range of variation of cp "will 
usually he much greater than that of 0, so that it appears reasonable 
to take Ap greater than A0. 

The boundary conditions in this region may appear in several forms. 
At a solid boundary, the condition is 


V 


-w 


d(p 


( 42 ) 


where 0-g = 0g(cp) is the equation of the boundary. On a line -where the 

region abuts a known hyperbolic region, u -will be specified. On a 
line of symmetry, the normal derivative of u -will be 0. 

The final, kind of boundary condition -will appear at a shock. The 
essential difficulty here is that the location of the shock is initially 
unknown. Ho-wever, the following method may be used to specify this con- 
dition. First, the position of the shock is estimated to as great accu- 
racy as possible. This may be done in a number of -ways. If the body is 
nearly the shape of another for -vdilch the shock is known (as a circular 
cone), the shock -will be close to -bhat known shock. If part of the 
shock is already kno-wn from calculations in an independent hyperbolic 
region, the remainder may probably be estimated fairly accurately. 
Finally, the method given by Li^thill in reference 21 may be used to 
approximate the position. In any case, the more accurately the shock 
position is known at the start, the easier will be the solution. 

Suppose now that the shock has been estimated in some manner. One 
then uses as the shock boundary condition the (known) values of the 
radial velocity u corresponding to that shock. Then the flow is 
calculated roughly and the -variation of u^, or approximately, v, is 
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obtained. At the shock this velocity must satisfy the known jun^) condi- 
tion in the normal velocity. From equations (31) to (33)^ this is 

Vg = -u sin ©(q + sin%) (43) 


Then, assuming 6 remains constant, the variation of Vg with 0 may 

he found (fig. 6). If one then plots also the variation of v from 

the calculated velocity field, extrapolating if necessary, a new shock 
position may he found. This relaxation of the shock position can 
readily he carried out concurrently with the relaxation of the field. 

Once this irrotational solution is found, the entropy dlstrihution 
at the shock is given hy 


r(r - l)s = logg 


\ 2yM. a±n9 cos% - (7 - l) 
7 + 1 


7 log. 


(7 + l)M^ sin% cos% 


2 2 2 

2 + (7 - l)M sin 0 cos B 


(i^4) 


This completes the description of the method of solution of the 
general conical- flow problem. It should he noted that, althou^ the 
relaxation procedure is applicable without having any initial idea of 
the solution, it will he difficult to solve unless a reasonable esti- 
mate of the field can he made at the start. 


EXAMPLE 


Consider a flat-plate, zero-thickness, triangular wing (fig. 7) 
a -12° angle of attack, having as the included angle of the leading 
edges, 90°^ and moving in an otheiwlse uniform stream at a Mach number 

of 3- Since the leeidlng edges are supersonic (the Mach angle is 

sin“^ 1/3 fti 20°), the top and bottom surfaces are independent. The 

lower surface is an expansion surface and has no shock and hence the 
flow on this side is Irrotational. The solution on this side will he 
discussed first (fig. 8). 
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T^om equation (24), 1/^2 ~ 2 .l 66 j and q = 16. 731°. Then, from 

tables of Prandtl-Meyer expansions, = 2.8902. Thus 


Vjj. = 0.8388 
, Vy = -0.0732 

V2 = 0.1783 



(45) 


These are the Cartesian velocity can^onenta in region CDF. Also, the 
equation of the start of the expansion BF is 


-cot d = 2.66 i 4 cos(cp - 11.740) + 0 . 9 ^n sin(cp - 11.740) (46) 

and the end of the expansion fan CF is given hy 

-cot e = 22.5608 cos(cp - 11.740) + 0.9577 Bin(cp - 11.74o) (47) 


The values for the Cartesian velocity components, given hy equa- 
tions (45), are constant in CDF. This is not the entire hyperbolic 
region. However, by using the def in ing equation of the characteristics 
(equation (I9)) and values of the velocity components calculated in the 
expansion fan, there was calculated that characteristic starting at B, 
which is the point of tangency of the Mach cone from the wing vertex 
and the start of the e3q)5Lnslon. This iq BCD. Althou^ there is a hyper- 
bolic region to the left of BCD, it is dependent on the unlmown elliptic 
one. Hence the remaining domain ABCDEA must be solved by relaxation. 

The boundary conditions are: 


ABCD ; u specified 


DE: 


^ 

sin 9 VVsin fi/ 




EA: 



(48) 


To get an Initial estimate of the solution, the linearized solution 
of the problem was found. The perturbation conq)onent of the radial 
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velocity 8u is given (in the elliptic region, that is, inside the 
Mach cone from the wing vertex) hy 


„ ua cos 0 

6u X 


-It 1 - tan 0 sin cp , -1 

— + — — — sin 


2 




1 - (m 2 _ l)tan 0 sin tp 


y 1 + tan 0 sin cp _i 
< + — sin ^ 

Vm2 - 2 


+ tan Q cos cp I tan" 


Vm^ - 1\/(1 - tan e sin (p)^ + tan% cos^(2 - 
1 + (m^ - l)tan 0 sin cp 


/m2 - 1^(1 + tan 0 sin cp)^ + tan20 cos2(p(2 - M^) 


2 tan G cos 


1 

cp)/l - (m^ - l)tan^0l 

-1 + tan20 

1 + (m 2 - l)cos\ 

J 


- It 


Then, the linearized solution is 


u = u cos 0 + 5 u 




The initial guesses for the values -of u were the values obtained 
from equation ( 49 ), modified in an arbitrary manner to fit the known 
boundary values and the actual shape of the elliptic (and dependent 
hyperbolic) region. 

The solution was then conq)leted by a relaxation process. Figure 8 
shows the configuration "together "with the various regions encountered. 
In figure 9 ^ "the surface pressure coefficient, given by 



is shown. For con^jarison purposes, the linearized pressure coefficient 
is also given. Because the ratio of angle of attack to Mach angle is 
hi^ (0.62), the linearized theory can hardly be expected to yield 
accurate results. This conpletes the solution on the expansion side of 
the wing. 
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On the con5)ression side, the solution in region PGH (fig. 8) ms 
found first. As before, = 2.166? and = 16.731°. Then, from 

the shock tables 1^2 = 1.52^. Thus 


= 0.7076 
Vy = 0.0607 

Vg = 0.15011- 






(51) 


The equation of the first part of the shock (PG) is then 
-cot e = -1.4318 cos(cp - 11.740) + 0.9577 sin(cp - U.74o) (52) 


As in the escpansion case, this analytic solution does not cover the 
entire hyperbolic region, but only that part limited by the character- 
istic GH, staxting from the point immediately downstream of the shock 

where = c^. As before, there is a small hyperbolic region 

which is dependent on the elliptic one. The remaining domain GHDEKG 
must be solved by relaocation. First, however,' the shock position is 
needed. The Mach angle is 19*5° ani the shock angle corresponding to a 
plane wedge at a 12° angle of attack and a Mach number of 3 is 29.5°. 
Hence the shock angle at the middle ms chosen arbitrarily as 27.5° aiii 
a smooth curve fitted in. Then, for the domain, these boundary condi- 
tions hold 


KGH: u specified 


HDE: Ug 




sin 0 


\lmf - 




(53) 


EK; 




= 0 


An Initial estimate of the solution ms made in the same manner as 
in the pre'vious case. Then, after the relaxation had been carried out 
so that the residuals were fairly small, the shock position ms recalcu- 
lated in the manner already described. This procedure was continued and . 
the solution obtained. Ho particular difficulties were encountered and, 
in this case at least, the shock position did not have to be changed a 
second time. 
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The entropy distribution along the shock was found to vary fairly 
slov^ly and hence no atteJi^jt was made to correct the solution for this 
variation. 

The pressure distribution on this surface is shown in figure 9 H^nd 
con^jared with the result predicted by linearized theory. 

In table I, the velocity field for the wing is presented. In the 
regions -vdiere the flow is plane, the Cartesian velocity con^onents are 
given and, in the remaining regions, the radial con5)onent is given. 


DISCUSSION 


The method described in this report for the solution of coniceQ. 
flows appears to be applicable to any such flow. With one possible 
exception, the relaxation procedure should converge fairly readily. That 
exception mlg^ht arise in the case where a large hyperbolic region where 
the flow properties varied rapidly is embedded in an elliptic region. 

An exan^jle of such a case is the tip region of a lifting wing with sub- 
sonic leading edges. In such a case there will be a very large upwash 
velocity at the tip (the linearized theory predicts infinite upwash). 

It also should be remembered that, as was found in the exan^jle 
discussed previously, the solution obtained for an independent hyper- 
bolic region will not be valid up to the parabolic line, but, rather, a 
portion of the hyperbolic domain will depend on the elliptic one. 

Regarding the calculated exan5>le, several things are noteworthy. 
First the pressure on the shock and expansion sides are both consider- 
ably hl^er than the linearized theory Indicates. Thus the lift of the 
wing is larger than that predicted by the linear theory. However, one 
would eipect that there will actually be separation on the expansion 
side of the airfoil and that therefore some of this added lift will not 
be realized. A second observation is that the extent of the constant- 
pressure region has been drastically changed. 


CONCLUSIONS 


A numerical method has been described for the calculation of the 
supersonic flow about cones for the case where the shock, if any, is 
attached to the cone so that the flow is conical. The procedure enploys 
a combination of the method of characteristics and the relaxation method. 
Also, an iterative process for obtaining the position of the shock is 
described. The method can be applied to any conical flow. The main 
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problem in practical application is the convergence of the relaxation 
solution. No difficulty on this score should he anticipated unless 
there is a lange embedded hyperbo3J.c region -where the flow velocities 
are expected to vary rapidly. 


Brown University 

Pro-vldence, R. I. , August 23, 1951 
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TABLE I 

VELOCITY FIELD IN DISTUEBED EEGION^ 


(a) In region PGH idiere = O.JO'jG, Vy. = O.O607, 
and = 0.1504. 


X 

0 

-10 

-20 

-30 

-40 

-50 

r60 

-70 

15 

18.5 
22 

25.5 
29 

32.5 

0.7438 

.7393 

.7318 

.7211 

0.7433 

.7388 

.7313 

.7208 

0.7418 

.7371 

.7296 

.7189 

0.7396 

.7344 

.7265 

.7158 





0. 7303 
.7219 
.7109 
.6975 

0.7250 

.7159 

.7043 

.6905 

.6752 





0.6958 

.6813 

.6654 







0.6545 







(b) In region CDF ^ere = O.8388, Ty = - 0 . 0732 , 
and Vg = 0.1783. 


X 

0 

-30 

-60 

-90 

-120 

-150 

-180 

0 

0.8278 

0.8278 

OJ 

CO 

• 

0.8278 

0.8278 

0.8278 

0,8278 

5 

.8365 

.8346 

.8299 

.8242 

.8189 

.8155 

.8144 

10 

.8405 

.8352 

.8252 

.8120 

.8029 

.7985 

.7973 

15 

20 



. 8 i 4 l 
— 

.7912 

.7809 

.7536 

.7778 

.7554 

.7774 

.7554 


(c) On characteristic BC.^ 


9 




-90 

0.8378 

-0.0682 

0.1607 

-105 

.8304 

-.0490 

.1000 

-120 

.8178 

-.0255 

.0470 

-135 

.8059 

-.0055 

.0100 

-148 

.8018 

0 

0 


^All velocities are radial unless othervd.se specified. 

%alues given are constant on rays extending throng 
wing tip. 
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(a) For shock attached to boTHidaiy surface tangent to body surface. 
PXane JACK is orthogonal to OG j plane JKCH is orthogonal to OK . 



OCj plane HKJ]^J2C is orthogonal to OK. 

Figure i;.- Constniction of solutions by method of characteid.stics for 
shock attached to body-tangent boundary surface and for esqjansion 
of flow in hyperbolic region. 
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FigTjre 7.- Flat-plate, zero- thickness, triangular wing with supersonic 
leading edges in uni f orm stream at Mach number of 3 - 



Figure 8 .- Construction of solution for flow about triangular wing with 

supersonic leading edges. 






Pressure coefficient 
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Figtire 9,- Stirface pressure coefficient. 
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